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We discuss the notion of quantum mechanical coherence in its connection with time evolution 
and stationarity. The transition from coherence to decoherence is examined in terms of an equation 
for the time dependence of the density matrix. It is explained how the decoherence rate parameter 
arising in this description is related to the "unitarity defect" in interactions with the environment 
as well as the growth in entropy of the system. Applications to the "Zeno- Watched Pot Effect" 
and gravitational interactions are given. Finally, some recent results on applications to macroscopic 
coherence with the rf SQUID, where the transition from quantum to classical behavior could be 
studied experimentally, are shown. 
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There's no doubt that the concept of time is absolutely fundamental, be it for philosophers, physicists or the man 
| and woman on the street. Kant even snubbed those haughty mathematicians and put time before number, making 
(S) ■ the intuition of time the origin of arithmetic [1]. Peculiar to physicists however, and probably not agreed upon or 
' well understood by these other demographic groups, is that with the 20th Century lessons of relativity and quantum 
C") , mechanics behind them, physicists in their overwhelming majority will agree that time should be viewed as a physical 
' process. There seems to be no avoiding-even in principle- the need to understand how the clock works as a physical 
23 ! object. 

We are familiar with the slowing down of time in special relativity, perhaps most easily understood as a transfor- 
i-£h I mation of fields due to velocity. The various fields change in such a way, according to relativity, that any physical 
process-including some "clock"- goes slower. Similarly, gravitational fields can affect the clock, as in the red shift. 
While these are interesting enough, there are further and perhaps less familiar aspects of time as a physical process 
■ resulting from quantum mechanics. This will be the subject of our discussion. 
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I. TIME EVOLUTION AS A COHERENCE PHENOMENON 



For a clock to work, according to quantum mechanics, it is necessary to have A) different energy states present, and 
B) coherence between these states. This might be easier to understand if we first look at the more familiar situation 
for space instead of time. In order to have a highly accurate definition of a point in space we need, according to 
quantum mechanics, a well localized wavepacket. This contains a large spread of momentum states. Indeed if we 
had only one momentum component the wavefunction would be completely spread out everywhere and in no definite 
position at all. So localization in position requires a spread of momentum states. 

But this is not enough. In addition these momentum states must be coherent. If one momentum state alone 
is spread out everywhere, it is only the phases between the different momentum states that allows them to cancel 
everywhere but in the vicinity of one point. And so if these states are added incoherently the system is still not 
localized. 

And so is it also with energy eigenstates and time. A single energy eigenstate by itself does not tic, it must "beat" 
or interfere with something else. Typically, with atoms or molecules, we have a second energy eigenstate nearby and 
we observe the oscillating electromagnetic field resulting from the beating of the two levels. Or we can suppose a 
i"T°-meson clock [2] or neutrino clock using the oscillations in flavor (strangeness, electron-muon number) resulting 
from the interference of different mass eigenstates. If the two atomic or molecular levels are not coherent, the emitted 
photons will be "chaotic" and the electromagnetic field has no well-defined phase suitable for a clock; if the two kaon 
or neutrino states [3] are incoherent there will be no oscillations. 

We thus see that in addition to the slowing down of time which is possible by relativistic effects, there is a further 
aspect to time as physical process: coherence and incoherence. 
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II. DENSITY MATRIX 



To deal with such questions of coherence and incoherence in quantum mechanics it is not enough to think about 
the wavefunction of the system in question. We need a quantity which can describe something like an average over 
different wavefunctions. 

The tool for this in quantum mechanics is the density matrix. The density matrix is easy to understand if we 
suppose having-in our minds-the complete, big wavefunction of the whole world, or at least that part of it that 
might interact with what we want to study. What we want to study will in general only be a subsystem of what is 
described by the big wavefunction and involves only a few of the total number of variables. Let us take the simple 
case where this subsystem is a two-state system, with basis states (j)i,(j) 2 . (One state would be too simple.) Then the 
big wavefunction can be written 

* = \x\,x 2 x N \ 1 > 4>! + \xi,x 2 x N ; 2 > (j> 2 , (1) 

where the x's are all the variables outside the subsystem, with N a big number like 10 23 , and the kets |.... > are the 
wavefunctions for these variables. The "1" or "2" indicate that these in general depend on the state of the subsystem. 
Now if we only look at the subsystem and not at the x's, that is are only concerned with operators that act on <j>i, tfi 2 , 
then all we need to know is the "reduced" density matrix of the two-state system. This is a 2 x 2 matrix which can 
be expressed [4] as follows 

_ / <xi,x 2 x N ;l\xi,x 2 x N ;l> <x\,x 2 Xn'A\xi,x 2 x N ;2> 

P \<xi,x 2 x N ]2\xi 1 x 2 x N ;l> <xi,x 2 Xn;2\xi,x 2 x N ;2> 

where the interpretation of this matrix can be given as: 

(how much of 1 coherence between 1 and 2 

(coherence between 1 and 2)* how much of 2 

where of course the "amounts" are probabilities and are normalized to one: ( how much of 1 + how much of 2 = 1). 
The justification for calling the off-diagonal elements the degree of coherence between the two states may be seen by 

considering two extreme cases: let \x±,x 2 xn\ 1 > and \xi,x 2 xjy;2 > be either the same or orthogonal. In the 

first case (f> can be factored out in Eq [1] and ^ is just a product wavefunction; </>iand (f> 2 are completely coherent. 
In the second case no operator working just on the (f>i, 4> 2 subspace can cause the two states to interfere, they are 
completely incoherent. If we think of the x's as describing an environment in which our subsystem is immersed, we 
can say that to the degree the subsystem is correlated with the environment, the subsystem is internally incoherent. 

Finally, this environment could be a "measuring apparatus". This is something that establishes a correlation- 
usually a strong one- between the state of the subsystem and the \x... >, which now represent the "apparatus". In 
the limit of very strong correlation we can know what 4> is just by looking at the \x... >, which is of course how 
a measurement instrument should function. On the other hand, we thus see how the measurement "collapses the 
wavefunction " of the subsystem: <px and (j) 2 become incoherent [4] when the |x... > states are very different. 





III. THERMAL STATE, STATIONARITY 

An easy and limiting case of incoherence to understand is the thermal state. The density matrix operator p for a 
thermal state is proportional to 

P ~ ^ HlkT , (4) 

where H is the hamiltonian operator, T the temperature and k the boltzmann constant. Looking at this expression you 
might be inclined to see only the usual boltzmann factor of statistical physics, which it is. However, more important 
for our discussion is the fact that it is an operator. Furthermore this operator is a function of the Hamiltonian operator 
only, T and k being just numerical parameters. This means that there is a definite basis in which p is diagonal, namely 
the basis of energy eigenstates \E >: 

Pe,e> ~ Ce5e,e> ■ (5) 

This diagonal property has the physical interpretation that different energy eigenstates are incoherent. That is, 
consider the expectation value of some operator O : 
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O = Tr[pO] ~ Ce6e,e>O e ,e> ~ C e O e ,e , 



(6) 



In words, the contributions to O for each energy state are to be added separately and incoherently. 

One obvious consequence of this incoherence between energy states is that nothing can be time dependent: as we well 
know the thermal state is stationary. That is, if we consider the time dependence of any operator 0{t) = e lHt Oe~ %Ht 
in the thermal state we get 

0(t) = Tr[p e im Oe~ im ] = Tr[pO] = constant , (7) 

since we can use the fact that e zHt commutes with p and the permutation property of the trace to bring the e lHt 
factors together and cancel them. 

Is the thermal state the only state like this? Obviously not, for Eq [7] we see that the important point is Eq [5], 
the diagonality, and not the precise value of the C E . The above arguments go through for any p = f(H) and not just 
Eq 4. For this reason we wrote a general coefficient Ce and not the boltzmann factor. So if it's not the thermal state, 
what is the characterizing feature of the density matrix [5]? The property which will be true of any state describe 
by Eq[5] is that it is stationary. Since there is no coherence between energy states no operator can show a time 
dependence. 

Conversely, if we know on physical grounds that some system is stationary, we know that its energy states are 
incoherent. An important case of this is the radiation coming from a thermal source. The radiation itself, say neutrinos 
from the sun, need not follow a boltzmann factor. But since the sun, (to a very high degree of approximation) is 
stationary, we know that different energy states are incoherent. Understanding this helped to clear up a long-standing 
confusion on the treatment of neutrino oscillations [5] . 



IV. LOSS OF COHERENCE 



Limiting cases like perfect coherence for a pure state, or incoherence of energy states for stationary systems, are 
easy to understand. However, we are often confronted with the more difficult situation of partial and time-dependent 
coherence. Typically we imagine the problem of a subsystem in a pure, perfectly coherent, state at t — 0; which as 
time goes on becomes more incoherent as it interacts with the external world. Although this is in principle a perfectly 
well defined problem: just go ahead, find the total wavefunction at any time and average out the unobserved variables 
by calculating Eq [2] — there is no simple way to do this perfectly rigorously for a many-variable system. We must 
resort to some kinds of plausible approximations or guesses. 

A simple and reasonable picture, from which however we can learn a lot, amounts to assuming that the outer world 
or environment has steady properties with time and there is a loss of coherence at a constant rate for our subsystem. 
Formally this is like the "golden rule" calculation used in decay problems or scattering theory. For the two state 
system as described by Eq [2] we can give a quite complete description of the phenomenology of the problem on this 
basis and draw a number of interesting conclusions. 



V. PARAMETERIZATION OF p 

To set up the framework for this approach, we first parameterize the 2x2 matrix for p in terms of the familiar pauli 
matrices a: 

P=i(/ + P.<7), (8) 

which is appropriately normalized to Trp = 1. The two-state subsystem is governed by an internal hamiltonian which 
completely describes its behavior in the absence of interaction with the environment. It is also a 2 x 2 matrix which 
we can parameterize in the same way: 

H = V-a , (9) 

where we have left out a possible constant term which simple represents an overall energy shift and has no interesting 
effects. 

These are just mathematical definitions, but because of the familiar identification of the pauli matrices with spin, 
they allow a simple intuitive picture. The 3- vector P is like a "polarization" . According to Eq [3] its vertical or 
"z" component represents the relative amounts of the two states "1" or "2". The "transverse" or x-y components, 



3 



according to Eq [3] represent the degree of coherence within the subsystem. In generalizing to more than two- state 
systems, this remains true; the diagonal elements represent the amount of the "quality" in question and the off- 
diagonal elements the coherence between them. Note that this is necessarily a basis-dependent picture: we must 
specify which "quality" we are talking about. 



VI. ISOLATED SYSTEM 



Naturally if the system is isolated and has no interactions with the outside world the situation is very simple. The 
evolution is governed by the above hamiltonian and the equation of motion for p given by the basic formula 

ip=[p,H] (10) 

Using the algebra of the a this translates into 

P=PxV. (11) 

Pictorially, the "polarization" P precesses around a "pseudo-magnetic field" V. Note that V can be time dependent 
if we wish. So far we haven't done very much, this is just equivalent to the evolution of a spin in a magnetic field. 
Observe that with this equation the length of P cannot change: 

dP 2 

— = 2P.(PxV) = 0. (12) 

That is, P rotates without shrinking. This represents the fact that a pure state remains a pure state for an isolated 
system. 



VII. ENVIRONMENT 



The non-trivial part comes in when we try to handle the interaction with the outside world. As mentioned, if we 
make the assumption of an outside world with constant properties, interacting repeatedly with our subsystem such 
that there is a steady loss of coherence, that is a constant rate of decrease of the off-diagonal elements, we can obtain 
a plausible generalization [6] of Eq [11]. This has the form 

P = PxV-flP T . (13) 

This involves one new parameter, D. This is the damping or decoherence rate and describes the rate of loss of 
coherence between the two basis states. P-r means the "transverse" or (x,y) components of P. If we mentally turn 
off the first term in the equation, the D term leads to an exponential decay of the overlap or coherence terms in p, as 
in the "golden rule", and this rate is independent of the time when we start. If our subsystem is immersed in a large 
essentially infinite system it is difficult to imagine any other behavior. 

We have made one inessential simplification in writing Eq [13]. In order to concentrate on the interesting interaction 
between the decoherence and the internal evolution we have left out a possible term D'P Z in which there would be a 
direct relaxation between the two basis states. For example in the problem of two states of a chiral molecule separated 
by a barrier [4] we are assuming that the temperature is low enough that there is no direct jumping over the barrier. 
This is the simple classical relaxation process which can easily be dealt with but is not essentially relevant to the 
quantum mechanical aspects on which we wish to focus. Similarly, it will be seen that at long times P in Eq [13] tends 
to zero. This may not be exactly correct if the internal hamiltonian is such that there is a constant energy difference 
between the two basis states, large compared to an ambient temperature. Then P z should end up being given by the 
boltzmann factor and not be zero. So we implicitly make the assumption that any constant component V z is small 
compared to the temperature. 

The most important aspect of the presence of the damping term is that now the length of P can shrink: 

dP 2 

— = 2P ■ P = 2P • P = -2DP ■ P T ^ , (14) 

in general, and so the density matrix can change its "purity". That is, even without explicitly "quality changing" 
relaxations the length of P will inevitably decrease, if there is a Pt component. Such a component however, will be 
produced, in general, by the "rotation" induced by V. 

Although we shall explain below how to arrive at a microscopic understanding of D, even the simple phenomeno- 
logical Eq [13] allows us to arrive at some interesting conclusions. One of these is the "Zeno Effect" . 
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VIII. ZENO- WATCHED POT EFFECT 



The Zeno or Watched Pot Effect, perhaps first proposed by Alan Turing [4], goes something like this. According 
to the "Copenhagen School"'s treatment of the "Measurement Problem" (all these things are in because I'm not 
sure they exist ), a measurement leaves the state of the system in one particular eigenstate of the measurement. For 
example after a position measurement, the object is in some definite position-not in a superposition of positions. So 
we or Turing might say, let's keep rapidly measuring the object's position. It will repeatedly be fixed to the spot and 
not be able to move at all! 

Sometimes this "paradox " is used to say there's something wrong with quantum mechanics or something of the 
sort. But actually it represents perfectly reasonable quantum mechanical behavior and in fact is just a solution to 
Eq [13], for the case of large D. 

To see this, consider the behavior of Eq [13] in the case of a constant V, say along the x-axis. If we start P 
along the z-axis, it will tend to rotate in a circle in z-y plane: P z -the probabilities or "amounts" we have defined as 
represented by P z will oscillate. Now if we turn on D, the oscillations will be damped and gradually die out. As we 
make D larger and larger the oscillations will become overdamped and disappear completely. Finally in the limit of 
very strong damping one finds [4] , [3] 

P z ~ e-#* . (15) 

In the limit of large D, P z is practically frozen and will hardly budge from the spot! According to this formula the 
time scale, initially set by the oscillation time l/V, gets stretched by a factor V/D, which can be enormous when it 
turns out D is something like the interaction rate in a gas (see formula for D below) and l/V is the tunneling time 
between two states of a chiral molecule. In this way we gave a second answer to Hund's paradox on the stability of 
the optical isomers [4] (the first was parity violation [7]). On a practical level this means that in relaxation processes 
at low temperature where the classical barrier-hopping type of mechanisms are frozen out and quantum tunneling 
takes over, we have "anti-intuitive" relaxation, where a bigger D means slower relaxation [4]. 

You might say the repeated "measurements" by the environment have fixed the state a la Zeno- Watched Pot but 
in any case it's just a result of a simple solution to Eq [13] and has nothing mysterious about it. 

At this point we should warn of a misunderstanding which sometimes arises. The arguments, or similar arguments 
leading to Eq [15] or its equivalent, depend very much on our dealing with a system with two, or in any case a finite 
number of, levels. They do not apply to a continuous density of levels, as for the decay of a quasi-stable state into 
the continuum. This misunderstanding has occasionally led some people to erroneously suggest that the decay of a 
particle or atomic or nuclear level could be inhibited by observing if it has decayed or not. This is obviously silly. 
Mathematically the difference between the two cases is that with two or finite number of levels P z has a "flat-top" 
near t = in the absence of damping. That is, P z <~ 1 — (Vt) 2 for small times; while for a true decay problem we 
have e~ rt ~ 1 — Ft, a linear behavior. This leads to qualitatively different behavior with respect to "stopping and 
restarting" the system as we do when we turn on the damping to get the Zeno- Watched Pot behavior. 

Another nice way of undcrtanding this (suggested to me by Michael Berry) is to consider the decay process as the 
tunneling through a barrier, as in the Gamow picture of alpha decay. Now, when there is tunneling through a barrier 
connecting two discrete states, as in the chiral molecule problem [4], the origin of the Zeno- Watched Pot effect may be 
understood as follows. Tunneling between two discrete states is a very delicate process and depends very much on the 
quasi-degeneracy of the energy between the two states. With E sp u t the energy splitting , the tunneling probability 
goes as ~ {utunnei/E split) 2 , with ^tunnel the tunneling energy. Since u tun nei is typically very small, so must E spHt 
be small to have a non-neglible tunneling probability. Now if some external influence is causing the energy levels 
to move around constantly, even by small amounts, say by some shifting of the two potential wells in the two-state 
tunneling problem, the degeneracy is constantly being lifted and the tunneling is slowed down. Depending on the 
ratio of D / 'u>tunnel it can be practically stopped. This leads to the Zeno- Watched Pot effect. But the situation is 
entirely different if we have a quasi-bound state decaying into a continuum of states, as in the alpha decay situation. 
Even if the energy of the quasi-bound state is fluctuating, it will always find some continuum state with which it is 
degenerate. Hence there is no Zeno- Watched Pot effect. 



IX. FORMULA FOR D OR THE "UNITARITY DEFICIT" 

We now come to the microscopic understanding of D, the damping or decoherence rate. Using our general ideas, 
we can derive a nice formula for this, which furnishes considerable insight into the "decoherence" or "measurement 
process" . 
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From its definition in Eq [13] and looking at Eq [3] we see that D is a rate parameter, one which gives the rate at 

which the overlap between |xi,X2 x^;2 > and \xi,X2 x^; 1 > is decreasing with time. We take our subsystem 

to be interacting with the big environment in a constant, steady manner-this is implicit in our assumption of a 
constant D with no further time dependence. We can model this situation by looking at the subsystem as being 
bombarded by a constant flux from the environment. This may be thought of as particles, or excitations like phonons 
or quasiparticles. Our basic idea is that if the two states 1 and 2 of the subsystem scatter this incoming flux in a way 
which is different for 1 and 2, then the overlap < x\, x-i Xn] l\x\,X2 Xn; 2 > decreases. 

We recall the concept of the 5 matrix of scattering theory. This is an operator which turns the incoming wavefunction 
into the outgoing wavefunction: 

4>out = Stpi„ . (16) 

Now the important point here is that states 1 and 2 may scatter the incoming wavefunction differently, so we have 
two 5 matrices, 5i and 52 [8]. If these arc different the incoming particle or excitation, which was uncorrelated with 
the state of the subsystem before the scattering, will be correlated with it afterwards. So every scattering decreases 
the overlap and if we work out what this means for D we get the following formula: 

D = (flux)lmi < i\(l - S{S 2 )\i > , (17) 

The expression is proportional to the incoming flux since it is a rate parameter and \i > refers to the incoming 
state, with a possible average implied if we have many kinds of incoming states. "Im " means imaginary part. There 
is also a significance to the real part, which is an energy shift of the subsystem induced by the environment [4]. 

Since usually we have S^S = 1 for the 5 matrix, the formula says that the decoherence rate is related to the 
"unitarity deficit" produced by the fact that the different components of the subsystem don't interact the same way 
with the outside world. 

Eq [17] has two illuminating limits: 

51 = 5*2 D = , no decoherence (18) 

If both states of the subsystem interact equally with the outer world, D = 0, there is no decoherence. This is in accord 
with our feelings about "measurement" . If the "outside" does not respond to the differences in the subsystem there 
is no "measurement" and no "reduction" of the state. However, we do not need to use this language, the equation 
stands by itself. Note an important point: interaction with the environment is possible without decoherence. 
The other interesting limit occurs if one state, say 1, doesn't interact so 5i = 1 then 

5i = 1 D = 1/2 x (scattering rate on state 2) (19) 

This result follows from an application of the optical theorem of scattering theory (cross section) <~ Im 5 and 
Rate = (flux) x (cross section). This corresponds to the naive expectation that the decoherence rate is the scattering 
rate. This is just one particular limit however, and there can be subtleties such as effects of phases on D [9]. 



X. THE UNITARITY DEFICIT AND ENTROPY PRODUCTION 



There ia a classical conundrum which says, in apparent contradiction to common sense, that the entropy of a 
isolated system cannot increase. This is particularly simple to show in quantum mechanics where the definition of 
entropy is — Tr[plnp\. Since for an isolated system the evolution of p is governed by a Hamiltonian, we get for the 
time dependence of the entropy, just as in Eq. (7), Tr\e lHt plnpe~ lHt ] — Tr[plnp\ = constant. This argument says 
that as long as the time dependence of p is governed by a Hamiltonian, the entropy is constant. 

This constancy would also apply to the entropy —Tr[p Inp] of our two-state system where p is given by Eq. (8) if 
it were isolated -not coupled to the environment. Or as we learn from the argument, it could even be coupled to the 
environment; but in such a way that the evolution of p is given by a Hamiltonian. However we see from Eq. (17) that 
the coupling to the environmnet is not governed by a single Hamiltonian but rather by two Hamiltonians, giving the 
two 5 matrices. If we had one Hamiltonian we would have 5i = 52, in which case D = 0, and there is no decoherence. 

Hence there is a connection between D and the entropy increase of the two state system. In fact diagonalizing 
Eq. (8) and taking the trace, we find for the entropy -Tr[plnp] = ln2 - |((1 + P)ln{\ + P) + (1 - P)ln(l - P)). 
Taking the time derivative, we find for the entropy change of the two state system 
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where we used Eq. (14). The s=s refers to the limit of small P. It seems intuitively clear that the rate of entropy 
increase and the decoherence rate should be closely related and Eq. (20) expresses this quantitatively. By appro- 
priate generalization of the various quantities, this could be extended to systems larger than just two states. Note 
P • Pt = P T is necessarily positive. Furthermore, in thermal equilibrium where there is no coherence between states 
of different energy i.e Pt — 0, there is no entropy production. 



XI. DECOHERENCE IN MESOSCOPIC DEVICES 



In recent years our subject has moved from the theoretical- philosophical to the domain of the almost practical 
with the realization of quantum behavior for large, essentially macroscopic devices- "mesoscopic systems". This has 
been given additional impetus in view of the possible use of such devices for the implementation of the "quantum 
computer" . Furthermore the subject is interesting in connection with the idea-to my mind wrong-that there might 
be some limit where large objects don't obey the rules of quantum mechanics. Decoherence is of course the main 
question for the observability of such effects and their possible practical use. 

One of the devices that has been studied in detail is the rf SQUID, where by suitable adjustment of parameters it 
is possible to bring the system into the configuration of the familiar double potential well separated by a tunneling 
barrier. The "x" coordinate stands for the flux in the superconducting ring, and the system obeys-in the absence of 
decoherence- a Schroedinger equation in this variable. The states localized in one of the potential wells represent the 
supercurrent flowing in a given sense around the ring, with a very large number (microamps) of electrons reversing 
direction when we go from one well to the other. Creation of quantum linear combinations of these states, which 
can occur by tunneling, would certainly be impressive evidence for the general applicabilty of quantum mechanics. 
Some beautiful experiments [10] [11] using microwave technique have seen evidence for such combinations in SQUID 
systems. 

We have suggested [12] a particularly simple way to both see the quantum linear combination of the two states and 
to measure the decoherence time of the system, hopefully directly sensitive only to the true quantum decoherence. 
It involves the idea of "adiabatic inversion" . This occurs when a slowly varying external field can cause a quantum 
system to reverse its state, as when a spin "follows" a rotating magnetic field and goes from up to down. This 
phenomenon is also often refered to as a "level crossing" . It is an intrinsically quantum mechanical phenomenon and, 
-important for us- is hindered when the decoherence time is short compared to the time in which the inversion takes 
place. 

We propose to produce such an inversion in the SQUID by sweeping an external field and then observing the 
reversal of the direction of the flux in the SQUID. Note that the system need not be "observed" until the procedure 
is over — our method is "noninvasive" . When the sweep is faster than the decoherence time the inversion should take 
place, and when it is slower it should be inhibited. We are witnessing the transition from quantum (tunneling allowed) 
to classical (tunneling forbidden) behavior as the decoherence is increased. 

Going from fast to slow sweeps, the sweep time where the inversion begins to become unsucccsful thus gives a 
determination of the decoherence time. A possible difficulty here is that the sweep cannot be too fast, otherwise 
the procedure becomes non-adiabatic. However, our estimates indicate that a region of SQUID parameters and 
temperature should exist where fast-enough sweeps are possible without violating adiabaticity. 

In order to study these points in more detail we have developed numerical simulations of such systems, both for 
the study of the adiabatic inversion (the logical NOT) as well as for a two-SQUID system operating as a CNOT 
quantum logic gate [13]. In Fig 1 we show the results of a simulation [14] for the one-SQUID adiabatic inversion. 
The decoherence time 1/D was chosen to be about 39 000 units and simulated as a random flux noise. The SQUID 
parameters were /3 = 1.19, L = 400pH,C = O.lpF, giving a time unit of 6.3 x 10 _12 s and so 1/D = 0.25/xs . This 
decoherence time would correspond to about T = 50mK in the estimate D = T/(e 2 R), with R — 5Mf2 [12]. The 
simulation included the first 8 quantum levels of the SQUID so that possible effects of non- adiabaticity are taken 
into account. 

The vertical axis in Fig 1 shows the probability for finding the flux in the ring reversed after the sweep. We 
see that while the inversion is essentially always successful for sweep times less than decoherence time, it becomes 
progressively less so for longer sweeps. Hence we are seeing the transition from quantum towards classical behavior, 
and a measuremnt of when this takes place furnishes a determination of the decoherence time. The gradual fall-off 
seen after the completion of the longest sweep is indicative of another difficulty, relaxation. Our final states will in 
general not be distributed according to thermal equilibrium, and the final detection should take place quickly on the 
time scale for relaxation. 
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XII. DECOHERENCE AND GRAVITY 



Although the above ideas are rather simple and phenomenological, they can be applied to a wide variety of interesting 
problems, where of course the various quantities like V and Si,^ must be adapted to each particular case. These 
range from understanding the permanence of optical or "crural" isomers (Hund's paradox) [4] , to the study of neutrino 
oscillations [3], to the mesoscopic devices [9], [12], and the design of quantum logic gates [13] just discussed. 

Here, however, I would like to conclude in a more speculative vein, coming back to the question of the need for 
coherence between energy states to define time. Is the loss of coherence always a more or less accidental happening, 
depending on the particulars of the situation? Or is there something universal about it? 

There is of course one kind of interaction which universally and inevitably couples to energy states and which in 
principle "measures" or distinguishes them: gravity [15]. Gravity couples to mass/energy and so states of different 
mass/energy interact differently. Indeed one can derive the gravitational redshift by considering the K° meson "clock" 
[2]. • 

For this clock, we have two mass/energy eigenstates and interference effects between them tell time. By the same 
token this must imply some degree of decoherence between the different mass/energy states due to gravitational 
interactions. There ought to be an imaginary or dispersive part to the redshift. Naturally because the coupling, 
Newton's constant, is so small we expect the effects to be negligible under everyday conditions. It is nevertheless 
amusing to see what the formulas look like. We shall incidentally find that there is in fact an "everyday" application, 
involving the passage of oscillating neutrinos near galaxies. 

Qualitatively, we expect in regions of large and rapidly varying gravitational fields that different energy states 
become incoherent. Since gravity couples universally, this cannot be dismissed as an incidental effect. It will affect 
all energy states, and thus time is slowed down equally for all processes. If the effects become strong, time in some 
sense becomes ill-defined or stands still. 

We try to model the decoherence due to gravity by calculating D for the K° clock. There are two interfering mass 
eigenstates and an environment interacting gravitationally with these two states. If we take this environment to be a 
flux of particles, we can use our formula Eq [17]. The calculation of the S matrices is (non relativistically) the same 
as for coulomb scattering, where the "charge" is the mass. The two components of the clock have different masses 
Mi, M 2 and so Si, 5*2 arc different and thus there is a contribution to D. In an impact parameter representation we 
have 

D = {flux) j 2nb db Im i(l - S\ {i (b)S M2 (b)) (21) 

with S(b) = e 2iS ^ and 

where a = GEM/v, with G the Newton constant and E the energy of the incoming particle. As is usual in coulomb 
problems there is a divergence at large distances, and so we introduce a large distance cutoff l m ax which must be 
interpreted according to the physical conditions. Taking the two S"s and the imaginary part we finally get 

D « (flux)l 2 max (Aa) 2 , (23) 

with Aa = G(AM)E/v, which we have taken to be small and expanded to lowest order. AM is (Mi — M 2 ). We can 
try to estimate this under thermal conditions, as in the early universe, where all dimensional quantities are related 
to the temperature T. We end up with 

D ^T 3 (GAM) 2 =T 3 ( ^ M f . (24) 

planck 

As would have been expected, T must be at least on the order of the planck scale for this to be significant. And then 
the decoherence rate is AM (AM / M p i anc k) , which is still small, unless we are considering AM also on the planck 
scale. The result Eq [24] is of order ^/M^ lanck since decoherence, as explained, is a unitarity effect and so comes 
in to second order in the coupling. Although conceptually interesting this seems very remote from any potentially 
observable effects. 

Nevertheless there is, as promised, a "present-day" application. We now know there are neutrino oscillations. So 
possibly if the time scale for the release of the neutrinos in the Big Bang were short compared to their oscillation time 
(then) , they would still be oscillating in step today. Instead of the K° clock we have a neutrino clock. A beautiful 
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clock-if we could ever read it-a precise timepiece going back to the very first minutes of the Big Bang. But there's 
a difficulty: in traveling to us the neutrinos will pass near various mass and gravitational field inhomogcncitics, like 
galaxies, and these will tend to mix up the phases. Because of the huge mass of a galaxy Aa will no longer be 
small. One finds [15], with typical mass parameters for neutrinos and galaxies, that essentially any neutrino traveling 
through or near a galaxy will be decohered. Unfortunately, there will be no oscillations to read. So much for the Big 
Bang Clock, but at least there's a "practical" application of decoherence theory involving gravity. 

Note that these ideas, based on the different interaction of different mass/energy states, do not contradict the 
classical equivalence principle according to which different masses move in the same way. The interaction is so cleverly 
constructed that while different masses follow the same path classically, they get different quantum phases [2], and 
these are essentially what concern us here. The interplay between the classical and quantum is subtle and interesting. 
The last has certainly not been said on these questions. 
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FIG. 1. Simulations of adiabatic inversion in the rf SQUID for sweep times long and short compared to the decoherence 
time. The vertical axis is the probability for finding the flux in the SQUID reversed. The simluation parameters have been 
chosen so that decoherence time is about 1 / D — 39 000 . For sweep times less than 1 / D the inversion is essentially completely 
succesful, while for longer sweep times it becomes less and less so . An experiment of this type would show the transition from 
quantum to classical behavior and permit a determination of the decoherence time. For the SQUID parameters given in the 
text the time interval corresponds to 6.3 x 10~ 12 s, so the decoherence time is about 0.25 us. 
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